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Abstract: 

In existing dark matter models with global symmetries the relic abundance of dark 
matter is either equal to that of anti-dark matter (thermal WIMP), or vastly larger, with 
essentially no remaining anti-dark matter (asymmetric dark matter). By exploring the 
consequences of a primordial asymmetry on the coupled dark matter and anti-dark matter 
Boltzmann equations we find large regions of parameter space that interpolate between 
these two extremes. Interestingly, this new asymmetric WIMP framework can accommo- 
date a wide range of dark matter masses and annihilation cross sections. The present-day 
dark matter population is typically asymmetric, but only weakly so, such that indirect 
signals of dark matter annihilation are not completely suppressed. We apply our results to 
existing models, noting that upcoming direct detection experiments will constrain a large 
region of the relevant parameter space. 
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1 Introduction 

Andrei Sakharov [1] pointed out that a primordial particle number asymmetry is produced 
given three ingredients: symmetry violation, C and CP violation, as well as a departure 
from thermal equilibrium. The observed baryon asymmetry suggests that these conditions 
have been met in the history of the Universe. Since global symmetry violation is expected 
to occur in the early Universe, we expect the existence of primordial asymmetries to be 
rather generic. 

The dominant contributions to the observed matter density in the Universe are due 
to dark matter and baryons. Their density fractions are remarkably well determined from 
cosmic microwave background data from WMAP7 [2]: 

n DM h 2 = 0.1109 ± 0.0056 n B h 2 = 0.02258 jftjjgj^. (1.1) 
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Yet, very little remains known about the nature of dark matter (DM). In most theoretical 
efforts, the problems of dark matter and baryogenesis are treated separately, with the 
implicit assumption that the comparable densities of these two types of matter is simply a 
coincidence. However, the fact that the abundances of baryons and DM are just a factor ~ 5 
apart may be an indication of a common underlying origin. The paradigm of asymmetric 
dark matter (ADM) has been suggested as a way of linking the asymmetries, and thus the 
abundances in the dark and visible sectors; see [3] and references therein . 

In ADM two distinct scenarios are typically considered [15-44]. In one case, a primor- 
dial asymmetry in one sector is transferred to the other sector. Here the primordial DM 
(vdm) arid baryonic (t]b) asymmetries generally satisfy either tjdm ~ Vb or i]dm "C r]B 
depending on whether the transfer mechanism decouples when the DM is relativistic or 
non-relativistic, respectively. Whereas if both asymmetries are generated by the same 
physical process then t)dm ^ Vb is also possible. Now, if, in analogy with the visible sec- 
tor, the present-day DM population is totally asymmetric, one can explain the coincidence 
expressed in (1.1) by tuning the DM mass to the value 

& DM VB M n 

m A DM = „ m p , (1.2) 

' Ib Vdm 

with m p denoting the proton mass. 

The final DM abundance is rendered asymmetric by removing the conventional sym- 
metric component that arises from thermal freeze-out. In the existing literature, this is 
typically achieved by the introduction of either a strong coupling, in analogy with QCD for 
the visible sector, or new light states. Both ingredients effectively induce an annihilation 
cross section sufficiently large to suppress the thermal symmetric component. Our results 
alleviate this hurdle to ADM model building and bolsters the case for their further study. 

Here we consider the coupled evolution of the symmetric and asymmetric populations 
via the Boltzmann equations and quantify how large the annihilation cross section needs to 
be in order for the present DM densities to be asymmetric. We will see that the symmetric 
population depends exponentially on the annihilation cross section, and hence that a weak 
scale force typically suffices to remove the symmetric component. In this sense, the weakly 
interacting massive particle (WIMP) thermal freeze-out scenario is not incompatible with 
the paradigm of ADM. 

We then discuss the implications for ADM and point out that the spectrum of possible 
ADM scenarios is much richer than previously thought. In particular, we claim that most 
weakly coupled extensions of the standard model are likely to interpolate between the 
extreme cases of purely asymmetric DM and symmetric DM. 

In this asymmetric WIMP scenario the present-day DM abundance is determined by 
a combination of its thermal annihilation cross section, its mass, and the primordial asym- 
metry. This relation among the short-distance and primordial parameters generalizes what 
is required in the extremely asymmetric or symmetric scenarios. In fact, for an asymmetric 
WIMP the total dark matter abundance is allowed to depend on a new observable, namely 
the ratio of the anti-dark matter to dark matter number densities. In contrast, in the 

1 There are also many other models seeking to link the abundances of baryons and DM [4-14]. 
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asymmetric scenarios considered so far this ratio is fixed to be zero, and consequently the 
present abundance is set by fixing the mass as in (1.2). Likewise, in the symmetric scenario 
the dark matter and anti-dark matter abundances are assumed to be exactly equal, and 
the present abundance is then set by fixing the annihilation cross section. For asymmetric 
WIMP dark matter, these two scenarios are recovered as limiting cases of either small or 
large present-day anti-DM particle abundances. 

We then apply our results to some existing models of ADM transfer operators [3], and 
investigate more generally the utility of weak-scale suppressed higher dimensional opera- 
tors. This is motivated by our findings that an asymmetric WIMP needs an annihilation 
cross-section only a few times larger than a picobarn to obtain the correct dark matter 
abundance. As an illustration, for the Higgs portal we find that requiring the correct 
annihilation cross-section leads to spin-independent direct detection rates at currently ob- 
servable levels. The next round of direct detection experiments can cover a substantial 
part of the parameter space. 

The outline of the paper is as follows. In Sec. 2 we solve the Boltzmann equations 
for a generic species in the presence of a primordial particle/anti-particle asymmetry. In 
Sec. 3 we specialize our results to the case in which the species is the DM, and discuss the 
implications of the asymmetric WIMP. In Sec. 4 we apply our results to existing scenarios 
of ADM and determine the limits imposed by direct detection experiments. We summarize 
our results in Sec. 5. 

In Appendix A we present accurate results for 2 — > n collision terms in the approxima- 
tion that the incoming particles are non-relativistic on average and the final state threshold 
is much larger than the masses of the incoming particles. In Appendix B we apply the 
results of the previous Appendix to the specific transfer operators considered in the text. 
Finally in Appendix C we compute the elastic scattering cross section for direct detection 
arising from the Higgs portal operator. 

2 On the origin of asymmetric species 

We begin with an analysis of the effect of a primordial particle-antiparticle asymmetry on 
a generic species X of mass m. Previous work on the relic abundance in the presence of 
an asymmetry appears in [45, 46]. 

We assume that the particle X is not self-conjugate 2 , and that a particle/anti-particle 
asymmetry in the A-number is generated at high temperatures. As the Universe expands, 
the number violating effects decouple at a temperature To and the asymmetry is frozen in 
for T < Tq. At this stage the number density of particles (n + ) and antiparticles (n~) is 
controlled by a set of coupled Boltzmann equations. Under reasonable assumptions 3 these 
reduce to: 

+ 311^ = -((Tannf) (n + rT - n^ q n^ q ) , (2.1) 

2 This is certainly the case if X carries a U(l)x global number, but our results apply more generally. 

s In writing (2.1) we assumed that i) our particle species is in a bath of particles in thermal equilibrium, 
ii) that no mass degeneracy between the two sectors is present, iii) that the dominant process changing the 
rv^ densities is annihilation, and iv) that the annihilation process occurs far from a resonant threshold [47] . 
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where (oarm^) = ao(T/m) n is the thermally averaged annihilation cross section times the 
relative particle velocity. In the absence of significant entropy production it is convenient 
to introduce the quantity = n^/s, where s denotes the total entropy density. In terms 
of this quantity, the particle-antiparticle asymmetry can be expressed by 

rj = Y+-Y~. (2.2) 

Without loss of generality, we conventionally define particles to be more abundant than 
anti-particles such that rj > 0. Notice that, as anticipated, r/ stays constant in the thermal 
evolution described by (2.1). 

The present densities turn out to be strong functions of the primordial asymmetry 
r\ and the annihilation cross section gq. We distinguish between two limiting scenarios: 

• In the first regime the dynamics can be considered "strong" and the final abundance 
is dominated by particles over antiparticles. In this case Y^ <C Y,£ ~ n. This is 
typical of most baryogenesis schemes if X is identified with the ordinary baryons. 

• In the second regime the dynamics is "weak" and the final abundances of particles 
and antiparticles are comparable, and separately much bigger than the asymmetry. 
In this case Y^ ~ Y^ ^> rj. This is generally assumed in the standard thermal WIMP 
scenario if X is identified with the dark matter. 

In the remainder of this section we will present a careful study of the Boltzmann 
equations (2.1) in the presence of a primordial asymmetry, and provide a quantitative 
assessment of what "strong" and "weak" dynamics mean (see Sec. 2.2). We will see that 
there is a broad range in parameter space between the two extremes discussed above in 
which the final abundances are comparable, namely 

Y--Y+- V . (2.3) 

2.1 The fractional asymmetry 

In order to solve the coupled system (2.1) it is useful to introduce the quantity 

rr Y~ , s 

Because the asymmetry rj is conserved, knowledge of the above quantity suffices to deter- 
mine the relative abundances at any time via 4 

Y + =n— — Y-=n-^—. (2.5) 
1 — r 1 — r 

Notice that by definition r satisfies < r < 1. 

We refer to r as the fractional asymmetry because it gives a measure of the effective 
particle-antiparticle asymmetry of a given species at any time. This should be compared to 



4 As anticipated, the following discussion does not depend on our convention that particles are more 
abundant that anti-particles. In particular, the physics is invariant under the transformation r/ — s> —rj, 
r — > l/r. 
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the primordial asymmetry rj, which provides such information only in the ultra-relativistic 
limit. The functional dependence of r on the primordial asymmetry rj and the strength 
of the particle interactions, {cr aQn v), will be given below. The parameter r is in principle 
an observable quantity as it controls the expected indirect signal coming from particle 
annihilation. Small r corresponds to an extremely asymmetric case in which indirect an- 
nihilation signals are small or even absent, whereas in the large r regime such signals are 
unsuppressed. 

We assume that the Universe is radiation dominated in the epoch of interest, in which 
case the Hubble parameter is 



Mpi 



with Mpi » 1.22 x 10 19 GeV. 

If the annihilation cross section for the species X is not too small compared to a 
typical weak scale process, the dynamical effects encoded in (2.1) become relevant when 
the particle is non-relativistic. In this case the equilibrium distributions can be taken to 
be 

/mT\ 3 / 2 
W ' 

In the above expression g counts the internal degrees of freedom and [i = £T is the chemical 
potential. The explicit expression of the entropy density per comoving volume 

2vr 2 



n 



eq 



n 



<■<! 



9 



-m/T 



(2.7) 



s(T) 



45 



(2.f 



allows us to introduce the convenient notation 

n, 



Y 



eq 



'eq 
S 



rn 



(2.9) 



with a = 45#/(4v / 2V/ 2 /i eff ). 

With these definitions the dynamical equation for r(x) following from (2.1) reads 



dr 
dx 



where 



and 



V45 



Y 2 

1 eq 



eq 



(1 



eq 



1/2 



1/2 



1/2 

9 C S 



Mpima , 

1 x d 9 c s \ 
4 g eS dx J ' 



(2.10) 



(2.11) 



(2.12) 



5 By requiring that sR 3 stays constant in time, taking into account the temperature dependence of Ii^q 
and using the relation t = t(T) given by the definition of H, one finds that dlog/igg- = fdlog g e Q- With 
the help of this result one verifies that (2.12) coincides with the expression used in [48]. 
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Figure 1. Evolution of Y (x) illustrating the effect of the asymmetry rj. After freeze-out both Y~ 
and Y + continue to evolve as the anti-particles find the particles and annihilate. The YsL curve 
shows the abundance for 77 = 0, a mass m = 10 GeV and annihilation cross-section ao = 2 pb. In 
contrast, with a non-zero asymmetry rj = t]b = 0.88 x 10~ 10 and same mass and cross-section, the 
more abundant species (here Y + ) is depleted less than when rj — 0. Also shown is the equilibrium 
solution Y eq (x). 

In the above we have also defined r eq = e~ 2 ^ x \ where £ is determined by 

2sinh£ = -^-. (2.13) 

* eq 

Notice from (2.12) that we have taken into account the temperature dependence in 
/i e ff and g e g. Because g e Q is monotonically increasing with T, we see that the parentheses 

1 12 

in the definition of gj is positive definite. In the numerical results that follow we use the 
data table from DarkSUSY [48] for the temperature dependence of g*(T) and h e Q(T). 6 

Eq. (2.10) reproduces the well known case 77 = for which one finds that r = 1 for 
any x. We will instead focus on scenarios with 77 7^ in the following. As shown in Fig. 1, 
the effect of nonzero 77 is to deplete the less abundant species more efficiently compared to 
rj = for the same annihilation cross section and mass. 

2.2 The relic abundance of asymmetric species 

Equation (2.10) can be solved by numerical methods and imposing an appropriate initial 
condition at a scale x = x% > 10, where the non-relativistic approximation works very well. 
Although we have chosen Xj = 10, we have checked that larger values (10 < X{ < Xf, where 
Xf is defined below) do not alter the final result. From (2.10) one sees that in the early 

6 Note that we use the notation of DarkSUSY for the massless degrees of freedom parameters. To 
translate our notation to that of Kolb and Turner [49] one should make the substitutions g e Q — > g* and 
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Universe r = r eq , provided the cross-section is not too small. In our numerical solutions 
r(xi) is chosen to equal its equilibrium value r(xj) = r eq (xi) = e~ 2 ^ Xi \ 

We now present a very accurate analytic approximation of the numerical results. The 
numerical, exact solution will be compared to it shortly. 

While there are no analytic solutions of (2.10), it is not difficult to guess the qualitative 
behavior of r{x). As already mentioned, for small x the solution is very close to the 
equilibrium solution r eq (x). As the Universe expands, r eventually has difficulty tracking 
the equilibrium expression which is exponentially decreasing. The easiest way to see this 
is to inspect the first equation in (2.10) and observe that the second term on the right-side 
is proportional to Y eq ~ e~ x . The deviation of r from the equilibrium solution begins to 
grow. 

The freeze-out scale x = xj is then defined by the condition that the two terms on 
the right hand side of (2.10) no longer balance each other, and dr/dx ~ dr eq /dx becomes 
comparable to the terms on the right-side of the Boltzmann equation. The behavior for 
x > xj is then dominated by the first term on the right hand side of (2.10) whatever the 
asymmetry is because of the exponential dependance of Y eq on x. We can therefore write 

r(x) = r eqJ e~ Xjl ^ m) (2.14) 

with 

$(x,m)= f dx' x'- n - 2 gV 2 (2.15) 

Jx f 

rm/xf dT f T y 



, 9* 

m/x m \ ni , 



and r eq j = r eq (xf). The present value of the fractional asymmetry is thus controlled by 
the quantities r eq j, Xf, m, and ultimately by £j. Importantly, observe from (2.14) that 
the fractional asymmetry is exponentially sensitive to the annihilation cross-section. 

It remains to find an estimate for the freeze-out scale. An analytic expression for Xf is 
obtained by the condition that dr/dx ~ dr eq /dx is comparable to either the first or second 
term in eq. (2.10). This leads to 

^^-SXrjgl^x— (2.16) 
or equivalently using the approximate relation (2.13), 



,tanh£ f = 1 ^ w S 3^, (2.17) 

with 5 a number to be determined by fitting the numerical solution. Because Xf turns out 
to be > 20 for the scenarios of interest (see the end of this section), we can ignore the 
3/2xf factor in the above expression. 
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As long as g]/^\rj < 2x 1 j +2 , the condition (2.17) can be solved iteratively. After one 
iteration, and with the help of equation (2.13), we find 

x f - log(^ ; a/ A) + - log " * J f , — 2 +■■■■ (2-18) 

2 log 2n+ WJja f \)-g*, f (5%) 

This expression is accurate at the percent level, so we will use it in our applications. Notice 
that the effect of a nonzero asymmetry rj on the freeze-out temperature is very mild (see 
for example Fig. 1). In fact, for all g^fXr] < 2x r j l+2 the scale Xf is very close to the result 
found in the symmetric limit, in which case 5 is typically taken to be 5 = n + 1 [46]. As 
already emphasized, however, the effect of rj on the anti-particle population is considerable. 
The asymmetry can only be considered small if g^fXr} <C 2x^ +2 , in which case one has 
~ 3> rj, as anticipated in the introduction. More generally, in the "weak" coupling 
regime glfXi] < 2x l J +2 we find that to a high accuracy r eq j ~ 1. The final estimate for 
the present fractional asymmetry in this limit thus reads 



if Sr < 1 

2x n + 2 



(2.19) 



where <3? is understood to be $ = <fr(oo,m). This result agrees with [45]. 

While it is clear that Eq. (2.19) describes the regime where < 1, it is less clear 
how far this approximation extends to the "strong" coupling regime in which r ro <C 1. In 
order to appreciate the range of validity of (2.19), we estimate the integral (2.15) assuming 

1/2 

that gj is approximately constant in the relevant range of temperatures. With this 
approximation one finds that $ ~ g l J"j /[(n + l)x^ +1 ], and (2.19) becomes 

I XV9 *J \ ( XT]9 *>f 2x f 1 ,oom 

' ^ " < XP 1 " (n + lK+i ) =6XP l-^p- (^TI) 1 • ,2 - U) 



To see how far the approximation extends, first note that for > O(10 -4 ), the exponent 
in (2.20) should be at most O(10). Since for a weak-scale cross section xj ~ 20, the upper 
bound on the exponent translates into an upper bound of Xr/g^ 2 /[2x^ +2 ] < 0(0.2), in 
agreement with our assumption (2.19). We thus conclude that (2.19) not only incorporates 
the intermediate limit anticipated in (2.3), but also the "strong" limit in which > 10 . 

As 2xj +2 — > Xrjg^J 2 the fractional asymmetry gets further suppressed since now r eq j 
can be much smaller than 1. Because in this regime r ro is essentially vanishing for all 
practical purposes, we will be mainly concerned with the regime (2.19) in the following 
sections. 

Finally, the estimate (2.19) turns out to be a very good approximation of the exact (nu- 
merical) result. We scanned a large range of parameter space ranging from r\ = O(10 ±4l )riB, 
m = O(10 ±4 )m p , and do several orders of magnitude above and below the "thermal WIMP" 
value cro t wiMP = 0(1) pb for both s-wave and p-wave processes. Based on this analysis, 
we have found that our analytic expression departs by at most 5% from the numerical 
result. Consistent with the derivation presented above, within the scanned range we have 
Xvtgl 1 ] /[2xf 2 ] < O(0.1). 
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3 Asymmetric WIMP dark matter 



In asymmetric dark matter (ADM) scenarios the dark matter is assumed to have a pri- 
mordial asymmetry rj. For definiteness we will measure the DM asymmetry in units of the 
baryon asymmetry via the relation 

V = er] B (3.1) 

where tjb = (0.88 ± 0.021) x 10~ 10 [50] 7 . The present-day abundance of baryons is ps = 
m p s t]b, whereas the DM abundance is 

PDM = m s(Y + + Y~) 

= ms(Y + -Y- + 2Y~) 

= m s ( rj + 2^-^- J , (3.2) 

where we have used (2.5) to express the present day density of anti-DM, Y~, in terms 
of rj and Tqo. The first term of (3.2) corresponds to the asymmetric component which 
survives in the limit — > 0, where the abundance scales as pdm ~ V- Whereas the 
second term corresponds to the symmetric component, in which case the abundance scales 
as pdm ~ l/icannv) when — > 1. In the intermediate regime, one may expect the total 
abundance to scale as simply the sum of these two approximate scalings. However, this is 
not the case in general since the symmetric component depends exponentially on 77 {cr ann v). 
We write the observational constraint (1.1) as 

m p V 1 + r oo J Mb 

with m p and m the proton and DM masses, respectively, and = Y^/Y^ the present 
fractional asymmetry. The present-day constraints (3.2) and ps = m p s rjs are equivalent 
to satisfying (3.3), generating the correct baryon asymmetry ijb, and having a theory for 
(3.1). 

The literature so far has focused on the following two extreme scenarios. In the first 
scenario (ADM) the DM is totally asymmetric, namely = 0. The observed relic abun- 
dance is then explained if the DM mass is tuned to the value 

atw ^ DM m P a A\ 

ADM: m ADM = S (3.4) 

\l B e 

and there is no strong constraint on the cross section. The cross section simply needs to 
be large enough (see below) to suppress r^. 

In the second extreme scenario (thermal WIMP) the DM is totally symmetric, namely 
Too = 1. From our analytic expression (2.19) we find that in this limit the relation (3.3) 
becomes 



Thermal WIMP: a , WIM p = 1 " B „ * J—. (3.5) 

r] B m p n DM Mpi^wjMP V vr 



' 1 T) and vb are normalized relative to the total entropy density (i.e., neutrinos + photons). 
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If <?* is assumed to be constant, <3? = gjj /{{n + ]> and the above formula becomes 

the one commonly used in the literature. The DM abundance is now determined by tuning 
the annihilation cross section, the dependence on m being very mild. 

For the intermediate scenarios in which < < 1 the physics interpolates between 
the above extremes. In Fig. 2 we plot the allowed parameter space in the m — <jq plane 
of a DM candidate satisfying the condition (3.3) for different choices of the asymmetry 
parameter e (see Eq. (3.1)). The two figures refer to purely s-wave and p-wave processes, 
respectively. For a given asymmetry e, we see that the allowed region is bounded from below 
by the WIMP cross section (Jq^wimp (3.5), and bounded from the right as a maximum on 
the mass at the ADM mass value (3.4). 

For e — > and m > 20 GeV the lines asymptote to a horizontal line representing the 
thermal WIMP scenario, where the cross section is fixed but the mass can vary. The regime 
m < 20 GeV is instead strongly sensitive to the abrupt change in the number of degrees 
of freedom around the QCD phase transition. The latter appears as a large gradient in 

1/2 

gj at T* < 1 GeV [48]. The effect on the relic abundance of the species is relevant 
if Tf = m/xf < T* , and can be seen in Fig. 2 as the bump at m < XfT* ~ 20 GeV. 
This change in the degrees of freedom also tends to reduce the exponent in (2.14) for light 
species. 

In Fig. 2 the fractional asymmetry varies from = in the upper part of the curves 
to — 1 when the curves overlap with the thermal WIMP scenario (e — 0). The same 
information is presented differently in Fig. 3, where we trade e for r^. In this figure the 
exponential sensitivity of to the annihilation cross-section is evident. 

A generic model for beyond the SM physics is expected to lie in between the two 
extreme scenarios described above. In fact as emphasized in the introduction, a primordial 
DM asymmetry (e 7^ 0) is generally expected; meanwhile most models beyond the SM 
naturally lead to WIMP-like cross sections. The resulting asymmetric WIMP scenario will 
typically have a > a G)W iMP and < 1. 

It is instructive to re-express the fractional asymmetry for an asymmetric DM species 
as a function of its thermally averaged cross section <jq and the quantity $ = 4>(oo,m) 
defined in (2.15): 



where uq^imp and ^wimp ~ 3> refer to a symmetric DM candidate of the same mass. 
The first equality in (3.6) is a good analytic approximation to the numerical solution 
(see (2.19)), whereas the second is exact and follows from eliminating r], A and m using 
(2.11) and the empirical formulas (3.3) and (3.5). 

From (3.6) we see explicitly that for WIMP-like cross-sections and larger, the depen- 
dence of Too on the cross section is very sensitive. For a more quantitative evaluation, we 
show in Fig. 4 the solution of the above implicit equation. Indeed, we see that an O(l) 
departure from the thermal WIMP cross section implies a significant change in r^. 

Indirect signals of annihilating asymmetric DM in the Universe are suppressed com- 




(3.6) 
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Figure 2. Here we plot the annihilation cross section cto required to reproduce the correct DM 
abundance floM via a s-wave process n — (above plot) and p-wave n = 1 (bottom plot) for a given 
dark matter mass m, and for various values of the primordial asymmetry r\ = er]B- The line for 
e = corresponds to the usual thermal WIMP scenario. Notice that the fractional asymmetry runs 
from Tea — in the upper part of the curves to = 1 when the lines converge on the standard 
thermal WIMP curve. The effect of the QCD phase transition appears as a bump at m < 20 
GeV, as anticipated in the text. Note that the bottom plot is basically enhanced by a factor 
^ n =o/^n=i ~ (n + l)xf compared to the former. As a reference, recall that 1 pb ~ 2.6 x 10~ 9 
GeV- 2 . 
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Figure 3. Same as in Figure 2 showing contours of constant r^, for s-wave process. The completely 
asymmetric scenario <C 1 corresponds to the top region of the plot having cross-sections always 
larger than the thermal WIMP cross-section. For reference we show the e = line which corresponds 
to the usual thermal WIMP scenario (roo = 0). 



pared to those of a symmetric thermal candidate of the same mass by the quantity 

xr^x (-?— V<1. (3.7) 



&Q,WIMP \i+r, 



One can understand the above formula as follows. The rate for indirect events is propor- 
tional to (cr ann v)r oo (nZ"/ ) 2 for an asymmetric candidate, or similarly to (<7 a nn , y)( r ^j = o) 2 
for a symmetric (thermal WIMP) candidate. Notice, however, that the number densities 
an< ^ n v=o are n °t t ne same. In fact, for DM candidates of the same mass it is the total 
(particle plus antiparticle) densities that must be equal, namely ri^, Q (l + r ro ) = 2n^ =0 . 
Taking the ratio of these rates gives (3.7). 

We plot the quantity (3.7) as a function of the cross section as a dashed line in Fig. 4. 
Interestingly, we find that potentially important indirect signals are still effective in the 
intermediate regime oq > (Jo,wimp m which r ro < 1. For example, a suppression of < 0.5 
is found for values of the fractional asymmetry below ~ 0.1 (i.e. for ao/cro,wiMP ^ 1-5), 
when the DM could already be considered asymmetric. 
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Figure 4. Present fractional asymmetry = Y^/Y^ (see (3.6)) for various values of the ther- 
mally averaged cross section <tq (defined as (cann^) = &o(T/m) n in the first section) in units of the 
value oq^wimp obtained for an exactly symmetric species of the same mass, see eq.(3.5). All the 
points in the curve, i.e. all solutions of (3.6), account for the observed DM density. The dashed 
line is the suppression factor (3.7) for indirect detection. 

4 Asymmetric WIMP dark matter: applications 

In the following subsections we apply our results to a few ADM models, choosing our 
examples from among those already existing in the literature [3, 36]. We demonstrate there 
is a large allowed region in these models in which the present-day fractional dark matter 
asymmetry is in the range < < 1 and, significantly, the dark matter annihilation 
cross-section in the early Universe can be comparable to the typical WIMP cross-section. 
Our results generalize the conclusions found in the previous literature, which only discusses 
completely asymmetric dark matter, corresponding to a present fractional asymmetry of 
Too = in our notation. 

In the examples we discuss the MSSM is assumed to couple to dark matter supernelds 
X, X via a heavy messenger sector. The dark matter chiral supernelds X, X carry ±1 charge 
under a hidden sector global U(l), are Standard Model gauge singlets, and have a SUSY 
invariant mass mj. For defmiteness we assume that after SUSY breaking the fermionic 
component tpx is the dark matter. Similar results hold for scalar DM candidates. 

The overall system is assumed to violate either lepton (L) number or baryon (B) 
number together with the dark (X) number, but to (perturbatively) preserve a linear 
combination of these symmetries. When the number-violating interactions are in thermal 
equilibrium, chemical equilibrium forces a relation between the asymmetries of the visible 
and dark sectors and therefore determine the quantity e introduced earlier (3.1). We will 
refer to this effect as a transfer mechanism. 
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4.1 Example 1: Lepton- number violating transfer operator 

In the first example we consider lepton number violation and take the leading transfer 
mechanism between the dark and leptonic sectors, after having integrated out the messen- 
gers, to be described by a nonrenormalizable term in the superpotential [3]: 



AW asym = XX ^ L , (1.1, 



with A taken to be close to the TeV scale. The operator (4.1) preserves a global L' = X—2L 
number, but violates both L and X. The quantum numbers L, X become good quantum 
numbers in the early Universe as soon as the transfer operator (4.1) freezes out, assumed to 
occur for simplicity at an instantaneous temperature T/> One can show that for A = 0(1) 
TeV Tq is below the scale T sp /j at which the sphalerons shut off, in which case the only 
source of lepton number violation at low temperatures is induced by (4.1). 

This model has either one or two stable particles, depending on the relative mass 
difference between the ipx^ and the lightest superpartner of the SM. To see that note that 
the model has a Z^ symmetry (see also [3]) which is the usual Z<i R— parity in the visible 
sector, plus a in the X sector that acts as +i on the fermionic component ipx and — i on 
the scalar component X (the X fields have the opposite charge). This symmetry remains 
unbroken by soft SUSY and electroweak symmetry breaking, and guarantees that at least 
one component of X is always stable. The other SUSY component of X generically decays 
through the transfer operator to the lightest component of X and SM fields. 

The transfer operator typically causes the lightest superpartner (LSP) of the SM to be 
unstable, if the channel is kinematically allowed. (The decay LSP — > XX + SM particles 
is consistent with the Z4 symmetry) . If m^sp < 2mx then this decay is forbidden and the 
model has two stable particles. In what follows we will assume misp > 2rax for simplicity. 

4.1.1 Chemical Potential Analysis 

To estimate the primordial DM asymmetry parameter e, see (3.1), we proceed in two steps. 
For a useful review on chemical potential analysis see [51]. 

As a first step we impose chemical equilibrium and charge neutrality at the scale T = 
T sp h > Td at which the sphalerons decouple, assumed to be below the Higgs condensate 
scale T c . At these temperatures the only dynamical particles are the DM and the Standard 
Model fields (including the top quark), and we have: 

Q oc 18/i u - I2fi d - 6/x e = (4.2) 

Mu ~~ Md = i^v — Me 
fi u + 2fl d + fl v = 
2/J,x + /J>v = 0, 

where the chemical potentials fj, u , fid, (J>v, Me> and fix refer to the Standard Model up, down 
quarks, the (purely left handed) neutrinos, the charged leptons, and the DM respectively. 
The first equation in (4.2) imposes electrical neutrality, the second accounts for the 
exchange, the third for the sphaleron process, and finally the last equation follows from 
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the transfer operator (4.1). The solution of the system (4.2) can be used to derive the 
primordial baryon and V asymmetries: 

VB = -yMe (4.3) 
25 n f(m/T sph ) \ 



= 1 y + 36 /(o) ; ^' 

where the function /(x) is defined by 

1 f°° v 2 

47r cosh 2 (jv^ + y 2 ) 

and accounts for the possibility that the DM be non-relativistic at the relevant tempera- 
ture [16]. If the DM is instead bosonic the cosh(x) function should be replaced by sinh(x). 

Below the sphaleron temperature the asymmetries for B and L' given above are con- 
served, while L and X are still violated by (4.1). In the second step of our chemical potential 
analysis we thus compute the present asymmetry for the X number at the temperature Td 
at which the transfer operator decouples. 

Chemical equilibrium at Tp imposes a number of conditions similar to (4.2). The two 
main differences are that at this scale the sphaleron process (see third line in (4.2)) is no 
longer effective, and that the top quark is not included. The system of equations for the 
chemical potentials is now solved by replacing the sphaleron process constraint with the 
initial conditions provided by (4.3), and not including the top quark contribution. The 
result is finally 

_ r]x _ / 309 + 22^g^ \ /(m/7b) 

~VB ^1026 + 72^^^1 /(0) ' { ' } 

Our result (4.5) reduces to the one found in [3] if the top quark is instead decoupled above 
the sphaleron scale, and if one assumes /(m/To) = f(m/T sp h) = /(0). We also checked 
that the result changes only mildly if one assumes that the massive vector bosons are 
decoupled at the scale Tjj. In all these cases one finds that (4.5) can be well approximated 
as e w 0.3/(m/2b)//(0). 

4.1.2 Transfer Operator Decoupling 

In order for the analysis of the previous sections to be applicable, the freeze-out temperature 
Tf must be below the scale To at which the transfer operator decouples 8 : only under this 
condition there exists a temperature range Tf < T < Td in which Eq. (2.1), and the 
subsequent analysis, is fully reliable. We will see in subsection 4.1.4 that, in typical UV 
completions of the model (4.1), this condition is ensured by higher dimensional operators 
that preserve the dark number and L, and which are therefore not involved in transferring 
asymmetries. 



8 We will comment on the limiting case Td = Tf in Section 4.1.3. For the moment, notice that the case 
Td < Tf is clearly meaningless: freeze-out always occurs when the strongest interaction decouples. 
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The Boltzmann equation governing the evolution of the X number at T > Td is given 
by Eq. (2.1), but with additional, A-number changing terms induced by (4.1). These latter 
involve the DM and the SM particles as final states, with typical reactions of the form 

The decoupling temperature Td is then estimated as the scale at which the strongest 
interaction mediated by the transfer operator - formally described by a "collision term" 
C on the right hand side of (2.1) - becomes comparable to the Hubble expansion rate, 
namely when 

C = H(T D )n x , (4.7) 

with nx the DM equilibrium number density. 

If the transfer operator decouples when the DM is relativistic, the Boltzmann suppres- 
sion in (4.5) is not effective and the parameter e does not depend on the DM mass. In this 
case the details of the transfer mechanism are irrelevant. 

The Boltzmann suppression in (4.5) is effective if the transfer operator decouples when 
the DM is non-relativistic. In this case only the precise value of Td, and hence the details 
of the transfer operator matter. In the specific case in which the neutralino is the LSP, 
we find that the dominant number- violating interaction mediated by (4.1) in the non- 
relativistic limit is associated to the process (4.6), and occurs via the production of an 
on-shell intermediate neutralino which eventually decays with a branching ratio BR = 1/2 
into an X- and L-number violating final state ipx^x ~~ ^x^x^- 

If we neglect the DM mass, then at leading order in T/m^ the collision term C reduces 
to (see Appendix B.l) 

C = feV ml ( — e-^' T . (4.8) 



4(2tt) 5 \mlAJ x v ,, x , 
In the above expression, g w is a weak gauge coupling, v u = v sin j3 is the vacuum expectation 
value of the up-type Higgs, and mo is the sneutrino mass. The exponential suppression 
in (4.8) reflects the fact that the process occurs far in the Boltzmann tail of the thermal 
distribution of the incoming DM particles. 

The factors of T in (4.8) arise from a partial integration in the process of taking the 
thermal average, and reflect the fact that the phase space integral vanishes when the energy 
of the DM particles is at threshold. Another way to understand these factors of T is to 
note that for the inverse process v\ —> ipx^x the powers of T arise solely from the number 
densities of the incoming particles simply because there is no threshold and the matrix 
element is set by the mass of the LSP. The inverse collision term is therefore proportional 
to nyn v ~ T 3 / 2 T 3 e" m */ T . 

The decoupling temperature is finally estimated using (4.7), with C given by (4.8) and 
nx the non-relativistic DM equilibrium density. We find that, up to negligible logarithmic 
corrections, Td is given by 

rn^-mx _ 1 ( g w v u V K ~ m x )M Pl m 7 f 

D ~ log7 ' 7 4(27r)3U?A; ^y/ 2 1/2 3/2 { -> 

\ 90 ) y e ff,D m x 
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LOD SB ™ 403 SB 

m DM [GcV] 

Figure 5. Contour plot for in the allowed parameter space for the model (4.1). Below the lower 
line and above the upper line the DM density is smaller and higher than experimentally observed, 
respectively. In the dark blue area 0.9 < < 1, in the blue area 0.1 < < 0.9, and in the light 
blue area < < 0.1. See the text for more details. Here we have taken rnp = A = 1 TeV and 
tan = 20. 

which approximately reads To ~ (m« — mx)/20 if the mass scales are of O(GeV-TeV). 
We now have all the ingredients required to estimate the primordial asymmetry parameter 
e, see eq. (4.5), as a function of the DM mass and the LSP mass mc. This estimate has 
been found here working under the assumption that the transfer operator decouples when 
the DM is already non-relativistic. Yet, because the Boltzmann suppression in (4.5) ceases 
to be effective when the non-relativistic limit is still a good approximation, it is clear that 
our estimate of e should also apply to the case in which the transfer operator decouples 
when the DM is relativistic. As emphasized above, indeed, in this latter case the details of 
the transfer mechanism do not matter. 

4.1.3 Discussion 

The main purpose of the present subsection is to scan the parameter space of the model (4.1) 
and show that there is a sizeable region in the allowed parameter space with a present-day 
fractional asymmetry in the range 0.1 < < 1. This is the asymmetric WIMP regime. 

The relevant parameters of the model (4.1) at low temperatures reduce to the decou- 
pling temperature Tp, the DM mass mx, and the annihilation cross section (encoded in 
Too). The requirement that the particle ipx saturates the observed DM density is imple- 
mented by (3.3) and it amounts to a constraint on these parameters. In Fig. 5 we show the 
resulting contour plot for the quantity defined by (3.3) as a function of mx and m^. 

The upper curve shows the line (3.4) at which r m = 0. In the area above this curve the 
DM abundance exceeds the observed value. This region of the parameter space is therefore 
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excluded. The lowest curve is determined by the condition To = Tf ~ mx/20, which 
approximately reads ~ 2mx in our case, at which r ro ~ 1 (and e <C 1). In the area 
below this latter curve the relic abundance of the species ipx is smaller than the observed 
value. This area is excluded here because ipx is assumed for simplicity to be the only 
DM candidate. The intermediate lines in Fig. 5 separate regions with different present-day 
fractional asymmetry r^. 

The plot can be understood as follows. For any r^, a given value for the DM mass 
completely determines e with (3.3). This fixes Tjj and in turn it determines a function 
m^(mx)- As long as the Boltzmann suppression in (4.5) is active the curve is approximately 
a straight line determined by mx/Tp ~ const > 1. This behavior is visible in the large 
DM mass regime of Fig. 5. The Boltzmann suppression shuts off when the previous lines 
hit the boundary of the non-relativistic regime mx/Tp ~ 1. For DM masses below this 
point the primordial asymmetry parameter is given by e ~ 0.3 and mx is fully determined 
by Tqo. This effect is seen as a vertical asymptote in the low DM mass regime of Fig. 5. 

It is evident from Figure 5 that a large region in parameter space is associated to the 
asymmetric WIMP scenario 0.1 < < 1. Here indirect signals are unsuppressed. 

4.1.4 Annihilation and direct detection 

In all the cases discussed above - either the strongly asymmetric C 1 or strongly 
symmetric ~ 1 limit - we see from Fig. 4 that the annihilation cross section can be 
of thermal WIMP magnitude. Such a relatively low cross section may be induced by 
annihilation into SM particles mediated by non-renormalizable operators suppressed by 
the very same scale A suppressing (4.1). From an effective field theory perspective in fact 
we expect that a generic UV completion of (4.1) also leads to direct couplings between the 
DM and the SM particles. For example, the following corrections to the superpotential 
and the Kahler are allowed by all the symmetries 

i2 XXH u H(i 
X 

AK sym = y 2 ^- 2 + ..., (4.11) 

and are therefore generic. 

We now show that both operators (4.10) and (4.11) provide annihilation cross sections 
of the right strength if A = O(TeV). These operators can therefore be responsible for 
keeping the DM in thermal equilibrium below Tp, and in particular for accounting for the 
values of r ro required to generate the correct relic abundance. As we will see, they can 
also provide direct detection signals. Let us now discuss a few annihilation modes in more 
detail. 

The correction to the Kahler triggers the process ipx'Px —> H, where I are SM leptons 
and ipx is a Dirac DM fermion. In order to avoid large flavor violating effects that might 
arise from (4.11), we follow [3] and assume that the only relevant coupling is with the third 
SM lepton generation. Neglecting the r and u T masses we find [3]: 

= <« 2 > 



AW sym = A 2 ^* + . . . (4.10) 
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It is easy to see that this channel can provide an annihilation cross section of the correct 
order for several DM masses. For example, for a light DM candidate, say with mx ~ 20 
GeV one should require A/y < 200 GeV, whereas for weak-scale DM masses, say mx ~ 200 
GeV one sees that A/y < 700 GeV suffices. If (4.10) provides the dominant annihilation 
mode, one finds no constraining direct detection bounds [3] . 

The dominant annihilation mode mediated by the operator XXH u Hd, on the other 
hand, occurs via the s-wave annihilation either to final state fermions - through the ex- 
change of the (virtual) CP odd scalar ^4° - or to a hA° final state. For the first channel 
the largest rate occurs if the final state fermions are top quarks, if kinematically allowed. 
Specifically, for ipx^x A * — > it and mx > mt, we have 

A 4 N c ml m\ \ m\ 

(fWj) = ta^S^r" A 2 {4m 2 x _ m 2 o)2 y 1 ~ ^T- ( 413 ) 

This can be of the correct magnitude to account for the DM abundance if ^4° is not too 
heavy, mx is not too small, and tan j3 not too large. For example, if we take mx = 200 GeV 
> m A o/2 and tan 2 /3 « 1 the averaged cross section is close to the WIMP value; for these 
values (a ann v) « A 4 x 4 x 10^ 9 GeV~ 2 (1.2 TeV/A) 2 . If this channel is kinematically closed 
the other channel may be open and provide a large enough rate. For 2mx > m^ + m A o, 
we have for the s-wave annihilation XX — > hA°, 



4 sin 2 (/3-a) / (m A o+m h ) 2 (m A o - m h ) 2 
{<T - V) = A 64vrA 2 V 1 ^nx~f ' ( } 



Unlike the previous cross-section, this one is not suppressed at large tan/3. If we take 
sin 2 (/3 — a) w 1, the averaged cross-section is (oa, im v) w A 4 x 4 x 10~ 9 GeV~ 2 (1.5 TeV/A) 2 . 

If (4.10) provides the dominant annihilation mode, the DM direct detection bounds are 
not generally negligible. In the absence of light degrees of freedom other than the MSSM 
particles and the DM, the "single nucleon" cross section for spin-independent elastic DM 
scattering mediated by virtual exchanges of CP-even Higgses is 

(lGeV\ 2 a N 

OXn = — -To , (4.15) 



where the spin-independent dark matter-nucleus cross section is 

a N = ^(Zf p + (A-Z)f n ) 2 (4.16) 

with f p and f n controlling the coupling between nucleons and the dark matter, and \it 
the DM-nucleus reduced mass. The nucleon coupling parameters f p and f n can be found 
in Appendix C and depend on the scale in Eq. (4.10) as well as parameters in the Higgs 
sector. 

Requiring that the annihilation cross section is larger than the WIMP value implies 
a lower bound on DM-nucleon cross section. Similarly, the bounds from direct detection 
can be translated into an upper bound for the annihilation cross section and thus a lower 
bound for r^. 
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Figure 6. Here we show the impact of the CDMS-II bound [52] axn = 5 x 10~ 44 cm 2 and the 
new XENON100 bounds [53] ax n = 3.3 x 10 -44 cm 2 on the parameter space of our model for 
m x = 400 GeV and A = 700 GeV (left plot) and A = 1000 GeV (right plot). To be consistent 
with the null results of CDMS-II, a parameter point must lie to the right of the red curve, while 
to be consistent with XENON100 bounds a point must lie to the right of the dashed blue curve. 
Obtaining the correct relic abundance requires (a ann v) > <Jq,wimp which occurs for points to the 
left of the solid black curve. For reference we include the dashed black curve corresponding to 
(cannw) = 2 (Jq^wimp, which is not present in the A = 1000 GeV plot. Note that the A — 1000 
GeV case is strongly constrained by direct detection bounds. 

To see that, let us assume that the dominant annihilation mode is XX — > hA° and that 
mx is at the weak scale. We also know from Fig. 2 that (dann^) ^ &o,wimp = 4.5 x 10 -9 
GeV -2 . The CDMS-II bound along with this constraint on the annihilation cross section is 
shown in Fig. 6. We also include the new XENON100 bound a Xn = 3.3 x 10 -43 cm 2 [53]. 
In Fig. 6 we have taken a dark matter mass mx = 400 GeV, a transfer operator scale 
A = 700 GeV (left plot) and A = 1000 GeV (right plot) and a light Higgs mass m h = 120 
GeV. We have also assumed m 2 A0 3> m 2 z such that the mixing angle a is determined at 
tree-level by a ~ f3 — ir/2 and the heavy Higgs mass is m\ ~ m 2 A0 . The relation between a 
and f3 implies that the annihilation cross section contours appear as straight vertical lines. 
The interested reader can find additional details in Appendix C. 

Our analysis of the annihilation modes has been model-independent so far. We end 
this section by emphasizing that, once a UV completion for the operator (4.1) is known, 
alternative reactions, and constraints, might arise. As an instance, consider the UV com- 
pletion proposed in [3]. One introduces two electro- weak doublets D,D with hypercharge 

and generalized lepton number L' = =pl, and two SM singlets X, X with L' = ±1. One 
readily sees that integrating out the D field one gets the operator (4.1) as well as (4.10) 
and (4.11). A typical feature of this UV completion is the occurrence of a mixing between 
the neutral component of the doublet D and the DM after electro- weak symmetry breaking. 
This mixing induces a coupling of the DM ipx to the Z° of order ~ 5A Uj rfU n ^/A, and thus 
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implies new annihilation and direct detection modes mediated by a virtual gauge boson. 
Given the lower bound on the annihilation cross section, one finds that the direct detection 
signal is always too large if the mixing-induced Z° exchange is the dominant channel. One 
should then suppress this channel by taking < O(10 _1 ). At this point the dominant 
annihilation mode would be given by (4.12). Alternatively one could consider the NMSSM, 
where a large contribution to (4.10) might be obtained integrating out the singlet S, also 
coupled to the DM via X'SXX. In this latter case the coupling of (4.10) would be naturally 
at the TeV scale, and freeze-out would thus be controlled by either ipx^x — > ^4°* — > it or 
XX — > hA°. Notice that in all of the cases discussed above the DM is assumed to be the 
lightest state beyond the MSSM. 

4.2 Example 2: Baryon-number violating transfer operator 

We now consider an example in which the transfer operator violates B and X, but not 
L. The model is again an extension of the MSSM, and it is described by the following 
nonrenormalizable term [3]: 



AW, 



asym 



XXu c d c d c 
A2 



(4.17) 



The superpotential operator (4.17) preserves a global B' = X + 2B number, but violates 
B and X separately. The physics of this model is very similar to the lepton violating 
operator (4.1), see for example the subsection 4.1. 

A chemical potential analysis similar to the one presented for the L-violating operator 
in this case gives: 



e = 



Vx 
VB 



1756 + 253 



f(m/T sph ) 
/(0) 



f{m/T t 



D, 



5472 + 836 



f{m/T sph ) 

TP 



- 50 



f{m/T D ) 
1(0) 



/(0) 



(4.18) 



Similarly to (4.5), this result is very robust and approximately reads e ~ 0.3/(m/To)//(0). 

The typical processes transferring the asymmetry between the visible and the dark sec- 
tors can occur through on-shell or off-shell neutralinos or squarks. The processes involving 
neutralinos are always strongly suppressed by phase space, and we find that on-shell pro- 
duction of squarks dominates for typical masses; see Appendix B.2 for details. Using again 
our approximate expression (A. 5) the collision term for the process ipxipx udd is found 
to be 



C 



9^vr 3 / 2 Ag 8 

Ill ~ 

32A 4 (2tt) 7 9 



T 



9/2 



1 + 



T m x 



ma rrin 



(4.19) 



where the powers of T have the same origin as the ones in (4.8), and we assumed that the 
squarks are degenerate for simplicity. 

The analysis can now proceed similarly to the lepton-number violating transfer opera- 
tor. The result can again be summarized with a plot analogous to Fig. 5, which for brevity 
we do not show. The comments in Section 4.1 regarding s-wave annihilation and direct 
detection via the superpotential operator XXH u Hd also apply in this case. 



- 21 - 



5 Conclusions 



In this paper we studied the implication of a nontrivial DM primordial asymmetry 77 on the 
DM relic abundance. For any 77, the physical quantity controlling the effective asymmetry of 
the species is the fractional asymmetry r(x) = Y~ /Y + . By solving the coupled Boltzmann 
equations for we found that r(x) depends exponentially on both the annihilation cross 
section and 77, see (2.19) [45]. The present anti-particle population of an asymmetric DM 
species is therefore extremely sensitive to the strength of the dynamics. 

The parameter space in a DM model with a primordial asymmetry 77 is defined as fol- 
lows. The DM mass m is bounded from above by the mass tiiadm of a totally asymmetric, 
Too = 0, DM species, see (3.4): 

m < m ADM = -77. m p . (5.1) 

For any 77 one finds m < vtladm as soon as one departs from the condition r^ = 0. In 
particular, very low masses m <C m p are always allowed provided — > 1. 

The thermally averaged annihilation cross section for an asymmetric DM candidate is 
instead bounded from below by the value found for a totally symmetric species of the same 
mass (thermal WIMP value), see (3.5): 

(crannf) > (o"ann« )wiMP- (5.2) 

For any rj the cross section approaches its lowest value when — > 1, while a larger 
cross section can still be accommodated in the regime r^ <c 1. Notice that because of 
the exponential dependence of r ro on the cross section, an effectively asymmetric DM 
population, say with < 10~( 2 ~ 3 ), is already obtained with an annihilation cross section 
of just a factor ~ (2 — 3) bigger than the thermal WIMP cross section. This latter point 
has interesting consequences for model-building, as it implies that new light states, often 
introduced in previous studies of asymmetric DM scenarios to induce large cross sections, 
are not necessary. 

A generic weakly coupled UV completion for the standard model will typically have 
weak scale annihilation cross sections. If the DM is stabilized by a symmetry distinguish- 
ing particles from anti-particles it is also likely that the DM has a nonzero primordial 
asymmetry 77. As a result, asymmetric WIMP scenarios are expected to be generic. These 
scenarios typically lie between the extremely asymmetric and symmetric limits. Surpris- 
ingly, the indirect signals from present-day annihilation can be much larger than what one 
would naively expect from an asymmetric scenario. 

To illustrate the generality of our conclusions we applied our results to two existing 
models of asymmetric DM [3]. We saw that the parameter space consistent with the 
present-day DM abundance and characterizing the asymmetric WIMP regime is sizable. 
This property stems from the fact that the present-day fractional asymmetry acts as 
new, low energy parameter. 

Moreover direct detection bounds impose interesting limits on the Higgs portal opera- 
tor XXHuHj. Although this operator was introduced to yield a sufficiently large annihi- 
lation cross section, the same interaction gives the dominant interaction with nuclei. We 
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found that the current bounds from CDMS-II and XENON100 constrain the model, but 
upcoming direct detection experiments will be sensitive to a wide range of the parameters 
relevant for the asymmetric WIMP through the Higgs portal. 
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A Collision terms with thresholds 

In this appendix we present a very general approximate expression for the collision term 
C for a generic 2- to A^-body process with an energy threshold M. For definiteness we 
focus on s-channel scattering, in which the cross section a (2 — > N) is a function of the 
Mandelstam variable s only. Our results can be straightforwardly generalized. 
In general, the collision term for a X\ + X 2 — > N process can be written as 

x n ' (/ wMe}) ( 2 *> 45(4, «, - EpJXIXI 2 ). 

where Ex 1 2 an d Px 1 2 are the energy and 4- momentum of the incoming particles, whereas 
Ef and p^ of the N final particles, (\A4 1 2 ) is the spin averaged matrix element, and fx lt2 are 
the momentum distribution functions of the Xj's. The step function enforces the condition 
that the CM invariant energy s of the incoming particles is above threshold, s > M 2 . 
After having integrated in the final state momenta, the above formula can be equivalently 
expressed as an integral of the cross section times the relative velocity a (2 — > N)v over the 
thermal distribution of the incoming particles. 

In our paper the energy threshold M is always a heavy particle mass which by assump- 
tion it satisfies M > 2mx (see the text). In this case one can see that the incoming DM 
particles must be highly energetic, px 1 2 — Ex x 2 ■ Hence we can write: 

C —I^I^'^"-^-^ (A - 2) 

/+1 roc poo E , E 

dcosO / dE 1 / dE 2 E l E 2 C{s/M 2 ) e~^^, 

-1 Jo Jm 2 /[2E 1 (1-cos8)] 

where s = 2E\E2{1 — cos 9) with 9 the angle between the incoming momenta, and the 
dimensionless function C = (E\E2v)a / (2 3 tt 4 ) = sa/ (27r) 4 has been defined for convenience. 
Notice that C is a Lorentz invariant function of s/M 2 . 

We are interested on the physics at temperatures at which the incoming particles are 
non-relativistic, namely T < mx 1 2 • Our basic assumption is therefore that the relevant 
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temperature T satisfies M/T S> 1. Physically, the fact that near threshold scattering 
dominates the collision term allows us to considerably simplify the above integral by making 
use of the following approximation 



dy F{y)e-y = F ( - n \y min )e~y- 



1 + 



1 



Umin 



(A.3) 



where J 7 ^ is the lowest nonvanishing derivative of J 7 , and y m i n S> 1. The theorem (A.3) 
will be used to approximate the integral in E2 and in cos 9, where in our case the expansion 
parameter is l/y m in = T/M. Here we only sketch a few intermediate steps. Using (A.3) 
the collision term becomes, up to subleading terms in T/M, 



C 



-M 



n+l 



- 



M J 

4 



-1 />oo 

m-l 



dx [2(1 - x)] 7 



dee n e T 



2(l-x)e / 



(A.4) 



rp \ n+3/2 

' I y J CW(1) 



+1 
-1 



dx [2(1 - x)]"/ 2 - 5 / 4 e ! v 1 



where C^ n \l) is the lowest nonvanishing derivative of C at s = M 2 , and where in the 
last approximate equality we have used a saddle point approximation and again neglected 
0(T/M) terms. In the last step we make use of (A.3) to evaluate the integral in x = cos 9. 
And finally we have: 



C 



4 nil 



M 2 



M J 



n+5/2 



cW(l)e"7 



1 + 



M M 



(A.5) 



Eq. (A.5) is the formula used in the text to estimate the collision terms for the lepton 
and baryon violating transfer operators. We checked that (A.5) is always very close to 
the exact, numerical result, which is consistently approached in the limit M/T S> 1. The 
dominant source of uncertainty in (A.5) comes from the 0(mx 1 S /M) corrections. 

In the general case the cross section can depend also on the Mandelstam variable t, 
in which case the function C will also depend on cos 9 explicitly, and the very last step in 
deriving (A.5) would change accordingly. 



B Collision terms for transfer operators 

B.l Lepton violating transfer operator: XXH U L 

In the early Universe the transfer operator (4.1) keeps ipx in chemical equilibrium with 
on-shell sneutrinos {y) through the processes 

Mx v* (B.l) 

Throughout we assume the LSP is a neutralino \. Then the sneutrino has a large decay 
rate v — > x v - Once H drops below the decay rate of the sneutrino, the processes (B.l) 
become Boltzmann suppressed and ineffectual. Below this scale, however, the V'x's remain 
in chemical equilibrium with the neutrinos through the processes 



4>xi>x < — > 



(B.2) 



mediated by an off-shell v. These are described by an effective Lagrangian obtained by 
integrating out the sneutrino, 

C eff ~ ^ ^(Mx^x) + h.c. (B.3) 

For my- > 2mx the LSP is unstable and decays to ipx^xv and ipx^x^- Still, the processes 

ipxipx <-> ipxtPx 17 " ( B - 4 ) 

remain in chemical equilibrium. These reactions are generated through two types of colli- 
sions, that differ on the kinematics of the initial state particles: 

• The collision is energetic enough to produce an on-shell v. This process is more 
Boltzmann suppressed compared to the following process: 

• The collision ipx^x is energetic enough to produce a real LSP (x) + v, which imme- 
diately decays back to tpx^x^- 

• The collision is below threshold and the neutralino is off-shell. 
The last item is described by the effective Lagrangian [3] 

£e//^§-r-(V^x) 2 ^. (B.5) 

With the assumption that ipx is non-relativistic at To we find this off-shell sneutrino + 
off-shell neutralino process to be subdominant. Let us now turn to approximating the 
collision term describing real neutralino production. 

The forward and inverse processes (B.4) contribute a "collision term" to the Boltzmann 
equation for nx 

d ^ + 3Hn x = C[f i ] (B.6) 
The "collision term" for the T-invariant processes is in general 

3 

2j ^ J (2tt) 3 2E Xi J l"27rf 2£ Vi ./ (2ir)«2£. v ., ./ !2tt y>2E x 
d 3 p Ul f d?p V2 



c ^ = _ ( iV(4) f d3px i f d3px i f d3px * f d3px * 



(2ir) 3 2E Ul J (2tt) 3 2E U2 
(2vr)V 4 )(^pf)(|^| 2 ) (f Xl fx 2 ~ fvjvjxjx) (B-7) 



We have assumed the occupation numbers satisfy /j < 1. The factor of 1/2 3 occurs for 
both the forward and inverse processes and is for integrating over one eighth of phase space 
since the two X's, two X's and two F's are identical particles. The factor of 4 is present 
because the process (B.4) changes the X particle number by four units. We note that 
the process ipx^Px ^ ipxipxvi' does not contribute to the collision term since it does not 
change the X particle number. 
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If these reaction rates are large enough in the early Universe these processes are in 
chemical equilibrium and the collision term vanishes, implying the condition 



2^x + ^ = (B.8) 

on the chemical potentials. As the Universe cools these processes eventually fall out of 
chemical equilibrium. When that occurs the asymmetries become frozen in. 

Next we turn to showing that C = C^ x ^ x ^ v where C^ x ^, x ^p is the collision term for 
the forward process ipx^x \P- This is expected since the neutralino does not carry X 
number. We will then apply our results from Appendix A to obtain an accurate expression 
for C^, x ^, x ^p. 

By assumed T— invariance (which is true at tree- level), the amplitude for ipxipx < — > 
XV i — > ip x ip x T7u factorizes as 

^x^x^x^x™ = M i>xi>xx^ 2 _ m 2 T im „ T „ M Mxxy ( B - 9 ) 

X XX 

where q is the momentum-transfer along the neutralino line, is the total decay rate of 
the neutralino, and ■M^ x1 p x ^ l/ is the matrix element for ipxipx < — > XiP- 

As noted above, the collision term for this matrix element is dominated by two regimes. 
In the first, q 2 <C m~, and describes the off-shell production of the neutralino. It is 
characterized by the effective Lagrangian (B.5) obtained by integrating out the neutralino. 

The other regime corresponds to real production of the neutralino. It simplifies in 
the narrow-width approximation, which is well-justified since the neutralino only decays 
through the transfer operator. Although this regime is characterized by a large energy in 
the initial state to reach the kinematic threshold (and leads to a Boltzmann suppression), 
the rate is enhanced in the narrow-width approximation and potentially important. 

To describe the collision term for the real production of x, we make use of the relation, 

2 



2 — m% + im^Tx 



S(q-m~) (B.IO) 



m x T x 



which is valid near threshold and in the narrow-width approximation. 
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Then for the forward reaction, the integral over the part of final state phase space 
dominated by the neutralino pole simplifies to 

r _L(A\ f d3pX i f d3pX 2 f f f f ^ S( 2 _ 2x ^ 

23 W i {2kY2E x J {2nf2E X2 TXlJX2 J {2nY2E v J (2tt) 4 ° W *Wx 



x(2vr) 4 5 (4) (px 1 + Px 2 -9-Ppi)(I-^Vx^x^iI 2 ) 
d 3 px 3 f d 3 px 4 f d 3 p v 



{2Kf2E x J (2itY2E x J (2^2E U2 
r) 4 5 {4) {pv 2 + PX3 + PX4 - q){\M~^ x ^ xV2 

\ ( A\ ( d3pX i f d3pX * f f f d L P ^l [ d4q S(n 2 m 2^ 27rT (X^^X^X^) 

2 [V J (2nY2E x J (2irY2E X2 Jx ' JX2 J {2^2E Vl J (2^°^ ^ ] r~ 
x(2tt) 4 5 ( - 4) ( PXi +p X2 - q - Pv 1 )(\M Mx ^ 1 \ 2 ) 
= 2 C Mx ^BR(x -> 4x$x") (B.ll) 



where C^ x ^ x ^p is the collision term for the forward process ii x ip x — > x v - Namely, 

r - = -f-\ o\ f d3pXl ( d3px2 f d3pw f d3p * 

fefexP " \2j [ Z) J (2ttY2E Xi J (2nY2E X2 J (2^)^ J {2^f2E^ 

{2^\Y. p "n M ^ x x V ?)fxJx 2 (B.13) 

i 

For this process the spin-averaged cross-section is 

(\M Mx5l u\ 2 ) = A 2 (2 PXl • PXa )(2pz ■ pv) (B.14) 

4 



\ 2 {s-2m 2 x ){s-m 2 ~) (B.15) 



where A = g w v u /(Am~) and s = Q 2 = (p Xl +p X2 ) 2 . Then integral over the two-body final 
state phase space is trivial, since the matrix element only depends on the center-of-mass 
energy. Using the standard formula 



/ 



d^> 2 (Q;p a ,p b ) = ^5^2 [(Q 2 -(^a + m h ) 2 )(Q 2 -(m a -m h ) 2 )] 1/2 (B.16) 



gives 



/ 



d^ 2 (Q ]W ^){\M Mx ^ 2 ){2^\Y, pl t) 



-^(s- m l) 2 (l-2 r f)0(s- m l) (B.,7) 



Now all but one of the angular integrals over the initial state momenta are trivial. 

The approximations we use are the following. Since T < m x < m^/2, the colliding 
Y>x's must be highly energetic to reach the kinematic threshold. Notice that this does 
not contradict the statement that the i/'x's are non-relativistic at the temperature Tp; 
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the previous observation is tantamount to saying that the DM particles that trigger the 
process of interest live far out on the Boltzmann tail. In the integrals we can therefore set 
pXi = and neglect terms of order m x /s and order T/mx- With these approximations 
the collision term simplifies to 

C Mxxv = ~J d cos 9 J dE i J dE 2 Ei E 2 e- {El+E2)/T 9{s-ml) (B.18) 
A 2 

with 

s = 2E 1 E 2 (l-cos6) (B.19) 
This expression for the collision term is precisely of the form (A. 2) with M = and 



2 9 vr 5 x \m\ 



Hence, from (A. 5) we find (n = 2) 

^ A 2« , T ,9/2 



C ^xv = - — J^[—) e~ m * /T - (B.21) 
Finally, with (B.ll) we get (4.8). 

B.2 Baryon-violating transfer operator: XXu c d c d c 
From the transfer operator (4.17) we consider the process 

->• 3?. (B.22) 

This occurs at tree- level through the dimension 6 operator (4.17), but its thermally averaged 
cross-section has a Boltzmann suppression factor that is sensitive to the mass difference 
between squarks and X. Specifically, for this process C/(nxH) ~ Exp[— (3mq — mx)/T], 
where we have taken equal squark masses. 

The spin-averaged matrix element for this process is 

(|M| 2 ) spin avg , = ^-{P 2 - 2m\), (B.23) 

with P = px 1 + Px 2 ■ The 3-body phase space can be decomposed in the product of two 
2-body phase space integrals using the identity 

d<S> 3 (P;p uP2 ,p 3 ) =d^ 2 (q;Pi,P2) d ^2{P;q,P3)(2iT) 3 dq 2 , q = Pl + p 2 . (B.24) 

The collision term can be written in the general form (A.2), with 

at | rP-rriq . . . 

C = 32(27r C )7A4 dq yjtf - AmyP 2 -(q + mtf^/P 2 - (q - m 9 ~) 2 . (B.25) 
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To make contact with the notation used in Section A we stress that P 2 = s and M = 3rriq. 

The integral expression C cannot be solved analytically. Fortunately, though, we just 
need its lowest nonvanishing n-th derivative (with respect to s/M 2 ) evaluated at threshold. 
One can verify that n = 2 so that the result (A. 5) in this case finally reads 



CxX3q 



9V2VV 2 N c l 
32A 4 (2tt) 7 



T 



m,- 



9/2 



-3m q /T 



rriq rrig J 



(B.26) 



There are additional processes implied by the operator XXu c d c d c involving the pro- 
duction of on-shell neutralinos. These are however less Boltzmann suppressed but receive 
a large additional phase space suppression. For rriq > we find that the process (B.22) 
is more important. In the main body of the text we focus on this latter regime. 



C Direct detection from the XXH u H d operator 

In this section we collect standard formulae [55, 56] relevant for the contribution of the 
operator 

A 2 



AWsym = —XXH u Hd 



(C.l) 



to direct detection. This occurs through the exchange of the two CP even Higgs bosons h 
and H. 

In 4-component notation one has the following interactions between the (Dirac fermion) 
dark matter and the physical neutral Higgs bosons, 

v\ 2 - 



C 



int 



2A 



ipx^x [cos(a + (3)h + sin(a + (3)H] 



Next recall the couplings of the physical neutral Higgs bosons to quarks: 



qqH/h v g - n p 
Integrating out the Higgses gives 



—p.uu [h cos a + H sin a] + 



v cos /3 



dd [—h sin a + H cos a] 



(C.2) 



(C.3) 



(C.4) 



where 



fu 



2A 



mi 



cos(a + j3) 



cos a 
sin/3 



+ 



H 



sin(a + (3) 



sin a 

sin f3 



A 2 

2A \ m 



-T7 [cos 2 (a) cot f3 — sin a cos a] H s- [sin 2 (a) cot j3 + sin a cos 

4 m H 



(C.5) 



A 2 ^ 
2A 

A 2 ^ 

2A V m 



mi 



cos(a + j3) 



sma 

cos (3 
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H 



sin(a + (3) 



cos a 
cos j3 



(C.6) 



[sin 2 (a) tan /3 — sin a cos a] H ^ [cos 2 (a) tan (3 + sin a cos a] ] 

^ m H J 
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Next we need the quark mass operator matrix elements 

rmA%> = (i\m q qq\i), i = p,n (C.7) 
From the literature one has [57]: 

proton : A$ ~ 0.023, A^ ] ~ 0.034, A® ~ 0.14, ^ t ~ 0.059 (C.8) 

neutron : A^ ~ 0.019, ~ 0.041, ~ 0.14, X^f t ^ 0.059 (C.9) 

Then the spin-independent target nucleus-dark-matter cross-section is 

a N = ^(Zf p + (A-Z)f n ) 2 (CIO) 

where \lt is the reduced mass for the target+dark matter system, and 

fi = m l Y J f q A<j\ i = p,n (C.ll) 
<? 

The "single nucleon" cross-section quoted in experimental limits is defined to be [58] 



0~Xn 



1 GeV\ 2 a N 



) J (C12) 



Note it is important to include the exchange of the heavy Higgs boson H since: i) its 
couplings to down- type quarks is enhanced by a factor of tan/3 at large tan/3; and ii) in 
the decoupling limit where the Goldstone bosons and light physical Higgs reside in H u , the 
heavy Higgs resides primarily in and its couplings to down-type quarks is proportional 
to cos a ~ 1. 
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